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ABSTRACT 



Mass fluxes J are computed for the extragalactic O stars investigated by Tramper et al. (2011; TSKK). For one early- 
type O star, computed and observed rates agree within errors. However, for two late-type O stars, theoretical mass-loss 
rates underpredict observed rates by ~ 1.6 dex, far exceeding observational errors. A likely cause of the discrepancy 
is overestimated observed rates due to the neglect of wind-clumping. A less likely but intriguing possibility is that, in 
observing O stars with Z/Zq ~ 1/7, TSKK have serendipitously discovered an additional mass-loss mechanism not 
evident in the spectra of Galactic O stars with powerful radiation-driven winds. Constraints on this unknown mechanism 
are discussed. 

In establishing that the discrepancies, if real, are inescapable for purely radiation-driven winds, failed searches for 
high-J solutions are reported and the importance of a numerical technique that cannot spuriously create or destroy 
momentum stressed. 

The Z-dependences of the computed rates for Z/Zq € (1/30, 2) show significant departures from a single power law, 
and these are attributed to curve-of-growth effects in the differentially expanding reversing layers. The best-fitting 
power-law exponents range from 0.68-0.97. 
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1. Introduction 

By observing O stars in nearby dwarf galaxies with metal- 
licities (Z) below that of the Small Magellanic Cloud, 
Tramper et al. (2011; TSKK) have extended the range of Z 
over which O-star winds have been detected and analysed. 
The resulting mass-loss rates ($) for Z/Zq ~ 1/7 are sur- 
prisingly high, leading the authors to question the Vink et 
al. (2001) scaling law $ oc Z a with a = 0.69 ± 0.10. 

Because of this result's importance for the later evolu- 
tionary phases of massive stars and perhaps also for the 
theory of stellar winds, models computed with the TSKK 
stars' effective temperatures (T e s) and surface gravities (g) 
are of interest in order to quantify discrepancies without 
relying on scaling laws. In addition, the validity of power- 
law scaling merits investigation. Accordingly, in this paper, 
mass fluxes J = Q/AirR 2 are computed for a wide range of 
Z at several points in (T g, <?)-space. 

Throughout this paper, the units of J and $ are gm s _1 
cm~ 2 and ©yr _1 , respectively. In these units, 



log® = log J + 2 log R/Rq - 3.015 



(1) 



2. Mass Fluxes 

In this section, the assumptions and procedures adopted 
when computing J from the theory of moving reversing lay- 
ers (RL's) are briefly reviewed, as are the curve-of-growth 
(CoG) effects that control J's dependence on Z . 



2.1. Assumptions 

The fundamental assumption is that the $ of a station- 
ary, radiatively-driven wind is determined by the constraint 
of regularity at the sonic point - i.e., that the transonic 
flow makes a smooth (analytic) transition from sub- to su- 
personic velocities. This assumption dates back to Lucy & 
Solomon's (1970; LS70) primitive theory of dynamic RL's 
and retains its role in updated versions (Lucy 2007, 2010a, b; 
L07, L10a,b). 

Additional assumptions are: 1) plane-parallel, isother- 
mal flow; 2) continuum emission occurs only at the lower 
boundary of the Schuster-Schwarzschild RL and has the 
frequency distribution of the TLUSTY model atmosphere 
with the prescribed T c g,g and Z (Lanz & Hubeny 2003) 
; 3) line formation occurs by non-coherent scattering with 
complete redistribution in the fluid frame; 4) line broaden- 
ing is due radiation damping and Doppler broadening with 
microturbulent velocity Vt — 10km s _1 ; 5) the RL's upper 
boundary condition of no incoming radiation is imposed at 
velocity v = 5a, where a is the isothermal sound speed. 



2.2. Procedures 

This multi-line transfer problem in a differentially expand- 
ing RL is treated with a Monte Carlo (MC) technique. This 
has the considerable merit of coding simplicity but has the 
drawback that gi, the MC estimate for the contribution 
of lines to the radiative acceleration, is not analytic. This 
poses the problem of how to determine the eigenvalue J 
from the regularity condition 
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g* - ge = at 



(2) 
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where = g — g e and g e is the radiative acceleration due 
to electron scattering. 

The procedure developed in L07 and LlOb replaces this 
point constraint by an integral constraint Q1.2 — 0, which 
is here conveniently rewritten as 



Al,: 



ge dx 



1 2 

r 



a 2 Inv + g*x 



(3) 



where the heights x\ and x-i correspond tou/a~l/2 and 
« 2, respectively. The interval (xi,x 2 ) includes the sonic 
point but excludes zones where the solution is affected by 
inexact boundary conditions. 

A MC simulation with AT indivisible photon packets is 
carried out in a model RL derived from the trial solution: 
gt = 9i at v i for z = 1,2, ... ,k, ... ,1 with gf. = at 
Vk = a. The resulting estimates ge then allow the left-hand 
side of Eq.(2) to be evaluated numerically, thereby testing if 
the equation is satisfied. With trial-and-error adjustments 
of gf for i ^ k, sequences of such simulations allow the 
eigenvalue J to be determined (L10a,b). 

This numerical procedure achieves consistency with the 
posed integro-differential problem in the limits / — > 00 and 
J\f — > 00, provided that the iterative adjustments of J and 
the gf have resulted in exact agreement between ge and gf 
throughout the RL. 

The main weakness at present of this technique is its re- 
liance on trial-and-error adjustment. In consequence, each 
solution demands a large commitment of the investigator's 
time. A reliable iterative algorithm is needed if large num- 
bers of solutions are required. 

Because the MC quanta are indivisible, none are cre- 
ated or destroyed within the RL. Accordingly, this MC 
technique does not spuriously create or destroy momen- 
tum. The sonic-point momentum flux Ja therefore derives 
from identifiable physical processes. 

2.3. Curve-of-growth effects 

Classically, the CoG refers to static RL's and relates the 
emergent line profile and its equivalent width to N, the 
absorbing ion's column density above the photosphere. 
However, the physical effects underlying the CoG are also 
relevant (LS70) to the J that can be radiatively-driven 
through the sonic point since ge at v — a depends on the 
co-moving radiation field, and this reflects line formation 
throughout the RL, including the quasi-static layers close 
to the photosphere. 

For simplicity, consider an outflow driven by one line 
from an atmosphere with continuum absorption due only 
to H and imagine varying Z with fixed T c g,g. We follow 
the discussion in LS70 for the A1548A component of the C 
IV doublet when the dominant ion is C III. For vanishingly 
small Z - i.e., negligible line formation, g\ oc n\, the ground 
state number density at = a. On the assumption of an 
ionization balance between radiative recombinations and 
photoionizations, n\ oc -n} c jn\ oc Z/p*. Accordingly, the 

sonic point constraint g\ = g* is obeyed at density oc Z, 
whence J = oc Z. 

In this weak- line limit, J oc Z because of increasing N . 
This effect is seen in Fig.2 of LS70 at logT cS « 4.55 and 
at w 4.25, where increasing N is due to shifting ionization 
balance as T e g increases. 



When increasing Z invalidates the weak-line limit, line 
formation reduces the co-moving flux seen by ions at v <~ a 
and so J(Z) flattens, reaching a maximum at J <~ J* = 
crT c 4 ff /c 2 , when all the photon momentum accessible with 
one line is converted into momentum of the transonic flow 
(LS70). In terms of the dimensionless rate <j) = ®c 2 /L, this 
maximum corresponds to <j> ~ 1. 

When increasing Z brings about the strong-line limit, 
corresponding to a line with a saturated Doppler core 
and well-developed damping wings, the photon momentum 
available at v ~ a is drastcally reduced, as therefore is J. 
This effect is seen in Fig.2 of LS70 at logT cS « 4.5 In this 
strong- line limit, the accessible photon momentum in the 
continuum is transferred to quasi-stationary gas (v 2 <C a 2 ) 
deep in the RL where it merely increases H p , the local scale 
height. 

With more than one line, the above CoG effects are 
superposed. However, in regions of high line-density, effects 
of overlapping line formation come into play. These defy 
simple interpretation and are the main reason for adopting 
MC methods. 



3. Models with varying Z 

In this section, models are computed at selected points in 
(T eS ,g, Z)-space. 



3.1. Parameters 

Because the RL's require input from the OSTAR2002 mod- 
els of Lanz & Hubcny (2003), their grid somewhat restricts 
the parameter choices made here. We adopt the TLUSTY 
labelling scheme, so that, for example, model Gt400g300 
has T c ff = 40.0kK, logg = 3 and Galactic metallicity. In 
order to test scaling laws over the full range of metal- 
licites found in the local Universe, a sequence of models 
with letter codes C,G,L,S,T and V are computed at each 
selected point (T e s,g). These correspond, respectively, to 
Z/Zq = 2, 1, 1/2, 1/5, 1/10 and 1/30, a range of 60 in 
metallicity. 

With regard to T eS and g, points in the OSTAR2002 
grid are chosen to match the best-fit parameters in Table 2 
of TSKK. The chosen points and the stars matched within 
errors are: 

t475g375 (IC 1613-A13) 

t425g375 

t350g350 (IC 1613-A15; IC 1613-C9; NGC 3109-20) 

t300g325 (WLM-A11; IC 1613-B11) 
The t425g375 sequence is included to avoid a large gap in 
T c ff. 

All models have v t = 10km s _1 . This canonical value is 
consistent with the highly uncertain TSKK estimates. Note 
that well-determined values of v t for Galactic O stars (e.g., 
Bourct et al. 2005) derive from UV spectra, which are not 
available for the TSKK stars. 



3.2. Model sequencies 

Values of J for the four sequencies are given in Table 1. 
The additional tabulated quantities are: 

Col 4- 4>i the mass-loss rate in units of L/c 2 . 
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Table 1. Computed mass fluxes J(gm s 1 cm 2 ). 



Seq. 




log J 


4> 




<(%) 




2 


-3.86 


429.8 


37.8 


80.0 




1 


-4.38 


129.8 


47.2 


81.5 


t475#375 


1/2 


-4.66 


68.1 


52.4 


80.3 




1/5 


-4.98 


32.6 


58.6 


69.9 




1/10 


-5.17 


21.1 


58.0 


58.7 




1/30 


-5.72 


5.9 


73.9 


23.6 




2 


-5.12 


36.9 


62.8 


89.4 




1 


-5.13 


36.0 


65.4 


86.4 


t425p375 


1/2 


-5.58 


12.8 


69.8 


64.6 




1/5 


-5.93 


5.7 


74.5 


38.7 




1/10 


-6.21 


3.0 


78.5 


21.4 




1/30 


-6.60 


1.2 


84.6 


9.9 




2 


-6.06 


9.2 


61.0 


64.6 




1 


-6.08 


8.8 


66.1 


53.3 


t350ff350 


1/2 


-6.10 


8.4 


70.0 


50.6 




1/5 


-6.78 


1.8 


76.4 


23.2 




i /in 
1/1U 


-0. 1 


1 Q 
l.O 


/».D 


14.0 




1/30 


-7.17 


0.7 


81.8 


5.7 




2 


-6.13 


14.5 


59.3 


50.4 




1 


-6.50 


6.2 


64.6 


33.2 


t300#325 


1/2 


-6.62 


4.7 


69.7 


21.8 




1/5 


-7.22 


1.2 


76.2 


5.2 




1/10 


-7.18 


1.3 


79.9 


4.1 




1/30 


-7.46 


0.7 


85.9 


2.1 
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Fig. 1. Mass fluxes J from Table 1 plotted against metallicity 
Z/Zq. The curves are labelled with the sequencies' T c g(kK). 

Table 2. Power law exponents a. 



Since <f> has a maximum value <~ 1 with one driving 
line (LS70), <p is the equivalent number of effective non- 
overlapping lines. 

Col. 5: rj, the percentage of the photon packets escaping 
across the RL's upper boundary that have not undergone 
a line scattering event within the RL. 

Since the MC quanta emitted at the lower boundary 
sample the continuum flux distribution of the TLUSTY 
model, a high value of rj indicates that much of this 
emission occurs in relatively line-free frequency intervals. 
Conversely, a low value indicates effective line-blocking of 
the star's continuum. 

Col. 6: (, the percentage of A^ 2 - see Eq.(3) - con- 
tributed by lines of Fe and Ni. 

This quantity measures the importance of Fe and Ni rel- 
ative to C,N,0,Ne,Si and S in accelerating matter through 
the sonic point. Note that the MC estimator for gi - Eq.(10) 
in L07 - readily allows the contributions of individual lines 
to be extracted. 

3.3. Dependence of J on metallicity 

In Fig.l, the J's given in Table 1 are plotted against Z/Zq. 
These log — log plots show significant departures from lin- 
earity and even include departures from the expected mono- 
tonic decline with decreasing Z. 

The quantities -q and ( provide a qualitative understand- 
ing of the J's dependence on Z. In the main, we see that J 
increases with Z but that 77 decreases. Thus, the enhanced 
J and (f> can be attributed to the shrinking of line-free fre- 
quency intervals as weak lines become significant absorbers. 



Seq. 


a ± a a 


CTlogJ 


£4750375 


0.97 ± 0.07 


0.10 


£425^375 


0.89 ± 0.07 


0.10 


£350^350 


0.68 ± 0.11 


0.17 


£300^325 


0.75 ± 0.10 


0.15 



The importance of CoG effects is illustrated by £• The 
spectra of Fe and Ni ions have numerous weak permit- 
ted lines arising from metastable levels. In contrast, the 
line-driving contributed by light and intermediate mass ele- 
ments is predominantly from strong ground-state resonance 
transitions. As Z increases, these ground-state transitions 
saturate, and their momentum tranfers then occur deep in 
the RL, increasing H p but contributing little to J. Because 
of this effect, the relative contribution of Fe and Ni to Ai ;2 
increases dramatically with Z. 

3.4. Power-law exponents a 

Conventionally, the Z-dependence of $ is approximated as 
a power law, so that 

log J = logJ G + a logZ/ZQ (4) 

Although Fig.l shows significant departures from this be- 
haviour, least squares fits of Eq.(4) to the data have been 
carried out and are given in Table 2. The intercepts logJc 
are not reported since they do not supersede the individual 
J-values for the G-models. 

In view of the significant departures from power-law be- 
haviour, the 3-D function J(T e f[,g; Z) will eventually have 
to be tabulated so that the evolution of massive stars can 
be followed as a function of initial metallicity. 
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4. Comparisons with observed &'s 

In this section, the $'s derived by TSKK for three of their 
stars are compared to the predictions of Sect. 3. The stars 
chosen are the one with high T e s and the two late-type O 
stars with the most accurate $'s. 

4.1. IC1613-A13 

For this star, TSKK estimate T cS (kK) = 47.6+4;^ and 

logg = 3.73+0 22 an d so tnc appropriate sequence from 
Table 1 is t475.g375. The metallicity of IC1613 is Z/Zq « 
1/7, which is bracketed by the models with Z/Zq = 1/5 
and 1/10. Logarithmic interpolation gives J = —5.07 
dex. Accordingly, with R/Rq = 11.4 from TSKK, $ = 
4ttR 2 J = —5.97 dex. This is to be compared with the 
TSKK estimate $ = -6.26toil dex. Thus theory over- 
predicts observation by 0.29 dex, but this is well within 
observational uncertainty and so there is no conflict. 

4.2. NGC 3109-20 

For this star, TSKK derive T eS (kK) = 34.2±|;^ and logg = 

3. 48+0^40 an d so the appropriate sequence is i350g350. The 
metallicity of NGC 3109 is Z/Zq w 1/7. Logarithmic inter- 
polation gives J = —6.77 dex. Accordingly, with R/Rq — 
24.7 from TSKK, $ = -7.00 dex. This is to be compared 
to the TSKK estimate that $ = — 5.41±g;|f dex. Thus the- 
ory undcrprcdicts observation by 1.59 dex, far beyond the 
observational uncertainty. 

In dimensionless terms, the TSKK rate is (f> — 76, in- 
dicating the number of effective driving lines required to 
accelerates the corresponding ,/ through the sonic point. 
Inspection of Table 1 shows that = 76 is only achieved 
for T c g ^ 45.0 and Z/Zq ^ 1/2, a domain excluded by 
TSKK's spectroscopic analysis. Moreover, the interpolated 
rj for this model is « 78 per cent, showing that there are ex- 
tensive frequency intervals not contributing to line driving 
of the transonic flow. 

4.3. WLM-A11 

For this star, TSKK derive T cS (kK) = 29.7+J;^ and lo 99 = 
3.25+0 19 an d so the appropriate sequence is t300<?325. 
The metallicity of WLM is again Z/Zq w 1/7. In this 
case, interpolation gives J = —7.20 dex. Accordingly, with 
R/Rq = 29.8 from TSKK, $ = -7.27 dex. This is to be 
compared with the TSKK estimate that $ = -5.56±ojo 
dex. Thus theory underpredicts observation by 1.71 dex, 
far outside the observational uncertainty. 

In dimensionless terms, the TSKK rate is <fi = 66 Thus, 
as for NGC 3109-20, this is only achieved for T cS £ 45.0 and 
Z/Zq ^ 1/2, a domain excluded by TSKK's spectroscopic 
analysis. The interpolated r\ for this model is also w 78 per 
cent, thus showing again that there are extensive frequency 
intervals that cannot be tapped to sustain a high J. 

5. Discussion 

In this section, the severe conflict between theory and ob- 
servation for the two late-type O stars is discussed. For 
clarity, the two extremes are considered. 



5.1. Observed rates in error 

The possible causes of overestimation of observed rates are 
discussed in detail by TSKK. As they fully recognize, a 
likely cause is wind clumping. Their estimates of <1> derive 

o 

from the partial filling-in of the He II A4686 A and Ha ab- 
sorption lines by wind emission - see their Fig. 2. Now, it 
has long been suspected that the instability of line-driven 
winds will lead to a clumpy outflow and hence enhanced 
emission in exactly these lines (Lucy 1975). Moreover, the 
standard diagnostic code CMFGEN (Hillier & Miller 1999) 
includes an adjustable clumping factor that demonstra- 
bly improves fits to observed spectra (e.g., Bourct et al. 
2008). Nevertheless, the model wind (Puis et al. 2005) that 
SDKK use to fit the data does not allow for dumpiness. 
Accordingly, one suspects that, if dumpiness were included, 
their genetic algorithm (Mokiem et al. 2005) would find a 
wide range of acceptable solutions: from high $'s with weak 
clumping to low $'s with strong clumping - i.e., degeneracy 
with respect to these two parameters. 

In view of the limited information content concerning 
winds in the optical spectra of O stars, the $'s of the TSKK 
stars are likely to remain highly uncertain until UV spec- 
tra can be obtained. However, in the meantime, as a control 
experiment, the TSKK methodology should be applied to 
optical spectra of Galactic O stars that already have anal- 
ysed UV spectra. If the $'s of Galactic O stars with similar 
spectral types to NGC 3109-20 and WLM-A11 can thus be 
derived with reasonable accuracy purely from their optical 
spectra, then the high TSKK $'s for these low-Z extra- 
galactic O stars would be more convincing. 

5.2. Theoretical rates in error 

The accuracy of the eigenvalues J is discussed in Sect. 4.2 
of LI 0b, from which it follows that errors ~ 1.6 dex would 
appear to be excluded. However, there is no proof that the 
eigenvalue J is unique. Accordingly, a search for a second, 
high-J solution has been carried out for St350g350 [Z = 
Zq/5). Thus, as described in Sect. 2.2, the vector g\ for 
i k is iteratively corrected with the aim of achieving 
dynamical consistency by matching the MC estimates gg. 
But now J is fixed at -5.19 dex, a value obtained by adding 
the discrepancy of 1.59 dex found in Sect. 4.2 for NGC 
3109-20 to the eigenvalue J = -6.78 dex from Table 1 (cf. 
Sect. 3.3 in LlOa). 

Fig. 2 illustrates the limited degree of consistency 
achieved when the search was abandoned. This may be 
compared with Fig. 2 in LlOb, where the successful solu- 
tion for Gt400g375 is plotted. In that case, J was not fixed 
but instead was iteratively corrected to obtain A^ 2 = 0. 

Fig. 2 shows that dynamical consistency is achievable 
for v/a Ss 0.6 but fails at higher velocities. Evidently, 
the constraint that the dynamical RL satisfies the regu- 
larity constraint (Eq. [2]) precludes adjustment into agree- 
ment with gg as the sonic point is approached. Accordingly, 
the inescapable conclusion is that a second solution with 
AlogJ ~ 1.6 does not exist. 

An identical experiment for St300g325 with J fixed at 
—7.22 + 1.71 = —5.11 dex also ended in a failure similar to 
Fig.2. 

The MC code cannot match these imposed sonic point 
momentum fluxes Ja because there are no physical pro- 
cesses available to provide these fluxes and the use of indi- 
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Fig. 2. Failed solution for 5**3503350 when J = -5.19 dex. The 
open circles are the gf for i — 1,2, ... ,k, ... ,1 from which which 
the model RL is derived. The filled circles are the MC estimates 
gi. The sonic point (i — k) is indicated. 



visible MC quanta precludes the spurious creation of mo- 
mentum (Sect. 2.2). 

5.3. An additional mass-loss mechanism? 

If the <& estimates of TSKK are accurate, then the above 
elimination of high- J solutions implies that an additional 
acceleration mechanism is operating besides radiative driv- 
ing. This is perhaps an unwelcome possibility, but where 
better to discover an additional mechanism than in the 
spectra O stars with low-Z's? 

For steady flow, Fig. 2 shows that this unknown mech- 
anism must provide an outward force per gm of ~ OAg at 
v = a in order to make up for g^s shortfall in satsfying 
Eq.(2). But extra driving is not needed at higher velocities. 
For v/a <; 1.5, Fig. 2 shows that a net outward acceleration 
is provided by §£. This is a consequence of matter being 
Doppler-shiftcd out of the absorption lines formed in the 
RL's quasi-static layers. 

The possibility that an extra mechanism may be re- 
quired only in the neighbourhood of the sonic point 
prompts the thought that small-scale, stochastic expul- 
sions may be accelerating blobs of photospheric matter to 
v/a <; 1.5, whereupon line-driving takes over and a high-$ 
supersonic wind ensues. 

A strong constraint on this unkown mechanism is pro- 
vided by the existence of a weak-wind domain for Galactic 
O stars (Marcolino et al. 2009; M09). From Table 3 in 
M09, we see that, observationally, this domain's extent is 
T eS {kK) e (31,34) and logg G (3.6,4.0). Thus, in (T cS ,g)- 
space, the two late-type O stars of TSKK are contiguous 
to the weak-wind domain, whose existence is solidly-based 
on UV spectra and is partially understood theoretically 
(L10a,b). Evidently, the unknown mechanism does not op- 
erate in the weak-wind domain. Accordingly, it must be 



narrowly focused in parameter space, perhaps also exclud- 
ing Z/Zq ~ 1. 



6. Conclusion 

The aim of this paper has been to test the theory of moving 
RL's at metallicities other than Galactic. The stimulus was 
the high empirical $'s reported recently by TSKK for ex- 
tragalactic O stars with Z/Zq « 1/7 and the doubts they 
raised about the Vink et al. (2001) scaling law. 

Although agreement within errors is found at 
T c ff(kK) ~ 47.5, discrepacies ~ 1.6 — 1.7 dex are found 
for two late type O stars with T e s(kK ) ~ 30 — 35. If these 
discrepancies are not due to wind clumping, an effect not in- 
cluded by TSKK, then an additional mass- loss mechanism 
is necessary, perhaps operating only at low velocities and in 
a restricted part of (T e g,g, Z)-space (Sect. 5.3). What this 
mechanism might be is an open question. 

Most probably, wind-clumping is the solution to these 
discrepancies, but the proof may not come until we can 
obtain spectra in the UV. 
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